PASJ: Publ. Astron. Soc. Japan , 1-??, — 
© 2010. Astronomical Society of Japan. 



Oh 



O 
O 



Study of the Intracluster and Intergalactic Medium 
in the Sculptor Supercluster with Suzaku 



(Received 2010 July 14; accepted 2010 September 15) 



Kosuke Sato, 1 Richard L. Kelley, 2 Yoh Takei, 3 Takayuki Tamura, 3 
Noriko Y. Yamasaki, 3 Takaya Ohashi, 4 Anjali Gupta, 5 and Massimiliano Galeazzi 5 
1 Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 
77 Massachusetts Ave., Cambridge, MA 02139, USA 
kjato@mit.edu 

2 NASAJGoddard Space Flight Center, Greenbelt, MD 20771, USA 
f - **) 3 Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency, 

£Nj ■ 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210 

4 Department of Physics, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397 
. 5 Physics Department, University of Miami, Coral Gables, FL 33124, USA 

in 

<N ; Abstract 

We studied the high temperature plasma in the direction of the Sculptor supercluster at z — 0.108 with Suzaku. 
' Suzaku carried out four observations in the supercluster: namely, A2811, A2811 offset, A2804, A2801 regions in 

Q^) \ 2005 Nov.-Dec, including the regions beyond the virial radii of these clusters. The study needed precise background 

estimation because the measured intensity of the redshifted lines, especially those from oxygen, were strongly af- 
fected by the the Galactic emission. The spectra taken in the regions outside of the virial radii of the member clusters 
were used as the background which included both the Galactic and Cosmic X-ray Background (CXB) components. 
We also used the background data which were taken near the Sculptor supercluster. Temperature and metal abun- 
(— > , dance profiles were determined to the virial radii of the member clusters, and then we searched for the oxygen line 

emission in the region outside of the virial radii of the clusters. As a result, the temperature of the clusters decreased 
toward the virial radii, and the spectral fits for the filament region did not require extra component other than the 
Galactic and CXB components. We constrained the intensities of Ovil and Ovm lines to be less than 8.1 and 5.1 
| xl0~ 8 photons cirT 2 s _1 arcmin~ 2 , respectively, as 2<r upper limits. The intensity of O VII indicates «h < 1.6 x 10~ 5 

OO 



crrr 3 (Z/0.1 Z Q r 1/2 (L/25 Mpc)~ 1/2 , which corresponds to an over density, 6 < 60 (Z/0.1 Z r 1/2 (L/25 Mpc)~ 1/2 . 
| Key words: galaxies: clusters: individual(A2801, A2804, and A281 1) - intergalactic medium - X-rays: diffuse 

Q\ . background 



1. Introduction ground emission and the intrinsic background. On the other 

hand, Chandra and XMM observations with grating spectrom- 

A large part of baryons in the local universe is thought to ex- eters also made the absorption-line study of the WHIM possi- 

ist along the large-scale structure as filaments of the intergalac- ble and restricted the WHIM density more significantly (Fang 

tic medium, called as warm-hot intergalactic medium (WHIM), et al. 2002; Fujimoto et al. 2004; Nicastro et al. 2005; Kaastra 

based on several N-body and gas dynamics simulations of cos- et al. 2006; Fang et al. 2007; Rasmussen et al. 2007; Takei et 

mological large-scale structure formation (e.g. Cen & Ostriker al. 2007a; Buote et al. 2009; Fang et al. 2010; Zappacosta et al. 

& \ 1999; Cen & Fang 2006; Cen & Ostriker 2006; Dave et al. 2010). 

2001; Chen et al. 2003). The WHIM, whose temperature and Because Suzaku XIS is characterized by lower background 

density is thought to be about T = 10 5 ~ 7 K and 10~ 6 -10~ 4 cm 4 , level and higher sensitivity below 1 keV (Koyama et al. 2007), 

is hard to detect with the present-day detectors such CCD cam- we can observe the intra-cluster medium (ICM) emission to 

eras. Thermal and chemical properties of the WHIM, if ob- the outer region of clusters near the virial radii (Sato et al. 

served, would bring us rich information about the structure for- 2007a; Sato et al. 2007b; Bautz et al. 2009; Reiprich et al. 

mation and evolution of the universe. 2009; Hoshino et al. 2010). Although Takei et al. (2007b) de- 

Recent Chandra and XMM-Newton observations with the tected no significant redshifted O lines from the outer region 

high resolution imaging spectroscopy enabled us to search of A22 18 with Suzaku, they set a tight constraint on the inten- 

for the WHIM emission (Soltan et al. 2002; Kaastra et al. sity (see also Tamura et al. 2008). Reliable estimation of the 

2003; Finoguenov et al. 2003; Galeazzi et al. 2007; Takei et foreground Galactic emission is most important in constrain- 

al. 2007a; Bregman 2007; Werner et al. 2008; Galeazzi et al. ing the WHIM emission. Past studies have shown that Suzaku 

2009). Werner et al. (2008) and Galeazzi et al. (2009) reported XIS can effectively constrain the Galactic emission (Gupta et 

the significant detection of the emission from the WHIM gas, al. 2009; Yoshino et al. 2009) better than the previous satellites, 

while the others could not because of the extremely low density The Sculptor supercluster (Schuecker & Ott 1991; Seitter 

of the gas compared with the high level of the the Galactic fore- 1992) at a redshift of z — 0.1 1 is one of the richest supercluster, 
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Table 1. Basic properties of the clusters. 
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redshift 


virial radius* 


Average temperature 


N ti 






(Mpc/arcmin) 


k{T) (keV) 


(cm -2 ) 


A2811 


0.1079 


1.96/16'5 


5.0 


1.54xl0 2U 


A2804 


0.108 


0.88/7'3 


1.0 


1.78xl0 20 


A2801 


0.108 


1.25/10'5 


2.0 


1.70xl0 20 



* At a redshift of z = 0.108, 1' corresponds to 119 kpc, and the virial radii, nso = 
1.95 h-^ a Vt<7">/10keV Mpc (Markevitch et al. 1998) 



Table 2. Suzaku Observation logs of the Sculptor supercluster. 



Region 


Seq. No. 




Obs. date 


(RA, Dec)* 


Exp. 


After screening 










J2000 


ksec 


(BI/FI) ksec 


A2811 


800005010 


2005- 


11-28T04:01:40 


(00 h 41 m 52. s 9, -28°33'19") 


30.6 


29.3/29.2 


A2811 offset 


800006010 


2005- 


11-28T21:10:16 


(00 h 40 m 42. s 7, -28°44'00") 


33.9 


32.4/32.3 


A2804 


800007010 


2005- 


12-02T08:51:05 


(00 h 39 m 35?2, -28°54'35") 


31.6 


30.0/29.8 


A2801 


800008010 


2005- 


11-27T08: 30:47 


(00 h 38 m 31?3, -29°02'48") 


34.8 


33.4/29.5 



* Average pointing direction of the XIS, written in the RA.NOM and DEC.NOM keywords of the event FITS files. 




h 42 m h 40 m h 38 m h 42 m h 40 m h 38' 



Fig. 1. Left: Combined XIS image of the four observations of the Sculptor supercluster in the 0.5-5.0 keV energy range. The observed XIS0, 
1, 2, and 3 images were added on the sky coordinate, and smoothed with tr = 16 pixel as 17" Gaussian. The instrumental background 
(NXB) were subtracted, and the exposure was corrected, though vignetting was not corrected. The green circles show the extracted region 
for the spectral analysis as shown in subsection 3.4. The subtracted point source regions are indicated by the white dashed circles. Right: 
Galaxy distribution associated with the Sculptor supercluster from 2dFGRS catalogue. The galaxies are selected between 0.103 < z < 0.113, 
and the nearer (z < 0.108) and farther (z > 0.108) galaxies than the sculptor redshift correspond to the solid circles and crosses, respec- 
tively. The blue boxes show the observation regions with Suzaku, and the green circle and dashed circles are same as the left panel. 



composed of more than 20 clusters of galaxies (Einasto et al. 
1997) including 5 Abell systems. It is located near the south 
galactic pole. Since WHIM is thought to distribute along the 
filamentary structures over several tens of Mpc scale connect- 
ing clusters of galaxies, it is likely that superclusters contain 
significant amount of WHIM. The Sculptor supercluster has 
been observed with ROSAT and ASCA. Spiekermann (1996) 
and Obayashi et al. (1998) found no indications of a largely ex- 
tended X-ray emission in their observed regions, which is over- 



lapping with our Suzaku field. On the other hand, Zappacosta 
et al. (2005) reported extended soft X-ray emission in 0. 1 -0.3 
keV band with ROSAT PSPC, and the region corresponds to 
~ 2-4 degrees east of our observations. They suggested that 
the diffuse emission was tracing the WHIM. 

Fang et al. (2010) reported a 3<x detection of an oxygen ab- 
sorption line from the WHIM in the "Sculptor Wall" at z ~ 0.03 
using the Chandra and XMM grating observations of a blazer 
behind this structure. The Sculptor Wall is located much 
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closer than the Sculptor supercluster in redshift and ~ 10 de- 
grees westward from our observations. The estimated hydro- 
gen density is n — 6 X 10~ 6 cirT 3 under the assumption of the 
uniform absorber distribution along the line of sight between 
0.028 < z < 0.032 using the O VII column density. Zappacosta 
et al. (2010) also reported results from the WHIM search in 
"Pisces-Cetus Supercluster" at z = 0.062 and "Farther Sculptor 
Wall" at z = 0.128 in the same line of sight as "Sculptor Wall". 
In contrast to the result of "Sculptor Wall", significant absorp- 
tions (cr > 3) were not detected from any of the strong H-like 
and He-like transitions of C, O, and Ne at the redshift. 

We use H = 70 km s -1 Mpc~' , f2 A = 1 - Q M = 0.73 in this pa- 
per. The redshifts of individual clusters are 0.1079 for A2811 
and the offset observations, 0.108 for A2801 and A2804, re- 
spectively, taken from NED as shown table 1. At a redshift 
of z = 0.108, 1' corresponds to 119 kpc, and the virial radii, 
r m = 1 .95 h-^ Vfc<r>/10keV Mpc (Markevitch et al. 1998), is 
1.96 Mpc (16.5') for an average temperature of k(T) = 5.0 keV 
of A2811, 0.88 Mpc (7.3') for k{T) = 1.0 keV of A2804, and 
1.25 Mpc (10.5') for k{T) = 2.0 keV of A2801. Because of this 
redshift, we can separate the WHIM O line from the superclus- 
ter relatively easily from the Galactic emission. Throughout 
this paper we adopt the Galactic hydrogen column density of 
N H = 1 .54, 1 .78 , 1 .70 x 1 20 crrT 2 (Dickey & Lockman 1 990) in 
the direction of A2811, A2801, A2804, respectively. Because 
the system is located near the south galactic pole, the Galactic 
absorption and emission are both low and have no complex 
feature. Unless noted otherwise, the solar abundance table is 
given by Anders & Grevesse (1989), and the errors are in the 
90% confidence region for a single interesting parameter. 

2. Observations and Data Reduction 

2.1. Observations 

Suzaku observed the Sculptor supercluster in November and 
December 2005. The observation logs are given in table 2, and 
the XIS image in 0.5-5.0 keV is shown in figure 1 left. Figure 1 
right shows a galaxy density map in 0.103 <z <0.1 13 using the 
2dF Galaxy Redshift Survey (2dFGRS) The XIS was oper- 
ated in the Normal clocking mode (8 s exposure per frame), 
with the standard 5x5 or 3 x 3 editing mode. During these 
observations, the significant effect of the Solar Wind Charge 
eXchange (SWCX) was not confirmed in ACE data 2 , while 
it was known that the SWCX affected some Suzaku spectra in 
the lower energy range as reported in Fujimoto et al. (2007), 
Gupta et al. (2009), Koutroumpa et al. (2009) and Yoshino et 
al. (2009). 

2.2. Data Reduction 

We used version 2.1 processing data, and the analysis was 
performed with HEAsoft version 6.6.3 and XSPEC 12.5.0ac. 
In the analysis of XIS data, we selected ELEVATION > 15° 
of the data set to remove stray-light from the day Earth limb. 
Event screening with cut-off rigidity (COR) was not performed 
in our data. Note that, even if we select the data with COR>8 
GV, the resultant fit parameters all stay within the statistical er- 

1 http://www.mso.anu.edu.au/2dFGRS/ 

2 http://www.srl.caltech.edu/ACE/ASC/ 



rors. In order to subtract the non X-ray background (NXB), we 
employed the dark Earth database by the "xisnxbgen" Ftools 
task. 

For the spectral fits of the ICM emission, we generated 
Ancillary Response Files (ARFs) for A2801/A2811 assumed 
the /3-model surface brightness profile, B = 0.8/0.82 and r c = 
0.'88/l.'6, respectively, reported in Obayashi et al. (1998), by 
the "xissimarfgen" Ftools task (Ishisaki et al. 2007). We also 
generated the ARFs for each observation assumed an uni- 
form sky of 20' x 20' for the Galactic and Cosmic X-ray 
Background emission. Note that we used the uniform ARFs 
for the spectral fits of A2804. We included the effect of the 
contaminations on the optical blocking filter of the XISs in the 
ARFs. Since the energy resolution also slowly degraded after 
the launch, due to radiation damage, this effect was included 
in the Redistribution Matrix File (RMF) by the "xisrmfgen" 
Ftools task. 

3. Results and Discussion 

3.1. Strategy of the spectral fits 

We analyzed the spectra of the Sculptor supercluster using 
the following way. 

1. Extraction of the point-like sources for each observation. 

2. Background estimation for the Galactic and Cosmic X- 
ray background (CXB). 

(a) Spectral fits of the outer region of the virial radii 
of the member clusters in A28 1 1 offset region as 
shown in figure 1 . 

(b) Estimation of the background models using the 
Suzaku observations near the Sculptor supercluster 
as shown in figure 3. 

(c) Comparing the above two results, we determine the 
background model. 

3. Spectral analysis of the ICM emission of A2801, A2804, 
A281 1 to the virial radius. 

4. Constraining the O emissions from WHIM in the fila- 
ment region between the clusters. 

3.2. Fore- and Background Estimation 

Using "wavdetect" tool in CIAO 3 , we searched for point-like 
sources in the Suzaku images in the energy range 0.5-2.0 keV. 
As shown by the white dashed circles in figure 1 , we detected 
and subtracted 7 point-like sources using 2 arcmin radius which 
corresponded to the HPD of the Suzaku XRT The flux level of 
the sources were > 1.5 x 10~ 12 erg cm~ 2 s _1 in 0.5-2.0 keV, and 
contribution from unresolved sources was taken into account 
when we subtract the CXB spectrum. After subtracting the 
point sources, we fitted the spectra integrated over the A28 1 1 
offset region excluding the virial region (r < r^o) of A2811, 
where a few galaxies existed as shown in figure 1 right, with 
the following models including the local hot bubble (LHB) and 
milky way halo (MWH) as the Galactic foregrounds and the 
CXB as a background: apecLHB+phabsx(apec MWH +pow CXB ). 
It was suggested that SWCX made a significant contribution to 

3 http://cxc.harvard.edu/ciao/ 
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Fig. 2. The panels show the resultant fits of the outer region of the virial radii of the member clusters with 
the models as shown in table 3. The data are plotted by black and red crosses for BI and FI, and the black 
and red lines show the best-fit model for the BI and FI spectra, respectively. The CXB components are 
shown in a magenta line, and the Galactic emissions are indicated by cyan and blue lines, respectively. 



Table 3. Summary of the parameters of the fits of the outer region of the virial radii of the clusters. We assumed zero redshift and a solar metallicity in 
the fits. 





case 1 


case 2 


case 3 


case 4 


case 5 


LHB 












kT (keV) 
Norm* (xlCT 3 ) 


10 +om 
z - uz -i.oo 


O.lO(fixed) 

2 61 +0 - 37 


<— (fixed) 

2 36 +0 ' 36 


<— (fixed) 

9 09+0.31 
Z - Jz -0.11 


<— (fixed) 

9 09+0.32 

Z - JZ -0.12 


MWH 












kT (keV) 
Norm* (xl(T 4 ) 


32 +0 03 


32 +0 03 
4 18 +8 - 67 


31 +o.03 
3 60 +5 - 63 


0.30(fixed) 

3 69 +05 ° 


<- (fixed) 

3 54 +0 M 


CXB 












r 

Norm*' + (xl(T 4 ) 


i 99+O.O8 
1,z -0.06 
7 q^+0.64 
'■^ -0.62 


i 99+O.O8 
1 - zz -0.08 
7 n 1+0.65 
'■^ -0.62 


1.4(fixed) 

o qn+0.48 


<— (fixed) 
9 01 +0 - 48 


<— (fixed) 
9.6 (fixed) 


X z /v 


218/167 


218/168 


231/169 


231/170 


235/171 



* Normalization of the apec component divided by the solid angle, SI", assumed in the uniform-sky ARF calcula- 
tion (20' radius), Norm = f n e n a dV /(4tt(1 + z) 2 D 2 )/f2 y xl(T 14 cm 40CJ7T arcmin , where Da is the angular 
distance to the source. 

' The referred CXB normalization in Kushino et al. (2002) is 9.6x10 for T = 1.4 in the above 
unit. 



Table 4. List of ^- 2 /dof in case 4 in table 3 for the fits of the nominal and considering the systematic errors such as the NXB level and contaminant of 
OBF. For details, see text. 



uncertainties 


nominal 


NXB 


contaminant 






+ 10% -10% 


+ 10% -10% 


*7dof 


231/170 


247/170 227/170 


231/170 231/170 



No. ] Suzaku observations of the Sculptor Supercluster 5 



Table 5. Suzaku Observation logs of NGC 253 offset and Q0109-3518. 



Region Seq. No. 


Obs. date 


(RA, Dec)* 
J2000 


Exp. 
ksec 


After screening 
(BI/FI) ksec 


Distance 1 
degree 


NGC 253 offset 803004010 
Q0109-3518 703037010 


2008-12-29 
2008-05-20 


(00 h 48 m 22. s 2,-25°02'51") 
(01 h ll m 45?0, -35°03'10") 


56.7 
30.0 


53.7/53.7 
26.0/25.9 


3.75 
9.01 



* Average pointing direction of the XIS, written in the RA.NOM and DEC.NOM keywords of the event FITS files. 
f Distance from A2811 




1234567 1234 



Fig. 3. Extracted regions of NGC 253 offset (Left) and Q0109-3518 (Right). Green boxes corresponds to the extracted region for the spectral fits, and 
the elliptical regions are eliminated from the spectral region. 



NGC253 offset region; case 2 NGC253 offset region; case 5 Q0109-3518 region; case 2 Q0109-3518 region; case 5 




Energy (keV) Energy (keV) Energy (keV) Energy (keV) 

Fig. 4. The spectra of NGC 253 offset and Q0109-3518. The notations in the figure are same as figure 2. 



Table 6. Summary of the parameters of the fits of the NGC 253 offset region and Q0109-3518 region. 



Region 


NGC 253 offset 




Q0109-3518 


N H (cm") 


1.42x10-"' 




1.93xl0 2U 




case 2 


case 5 


case 2 


case 5 


LHB 










kT (keV) 
Norm* (xl0~ 3 ) 


O.lO(fixed) 

9 4S+0-74 

^• HJ -0.78 


<- (fixed) 

2 56 +0 - 65 


<— (fixed) 

o qo+0.76 


<— (fixed) 

9 C7+0.72 
^■-"-0.72 


MWH 










kT (keV) 


n o 9 +0.08 
-Q.Q6 

7 CQ + I.I4 


0.3 (fixed) 


34 +0 07 

-0.06 


0.3 (fixed) 


Norm* (xl0~ 4 ) 


3 00 +0 - 72 




9 7Q+0.84 


CXB 










r 


i oo+0.09 

L - JO -om 


1.4 (fixed) 


1 30 +013 


1.4 (fixed) 


Norm*' + (Xl0~ 4 ) 


10 26 +0 <54 


9.6 (fixed) 


8 49+ 104 


9.6 (fixed) 


x z /y 


180/155 


185/158 


204/155 


207/158 



* Normalization of the apec component divided by the solid angle, tl l , assumed in the uniform-sky ARF calcula- 
tion (20' radius), Norm = f n e n n dV I (An{\ +i) 2 Dj)IQ^ xl(T 14 cm 400yT arcmin , where D\ is the angular 
distance to the source. 

' The referred CXB normalization in Kushino et al. (2002) is 9.6xl0~ 4 for T = 1.4 in the above 
unit. 
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unabsorbed emission-line component (LHB) even in the case 
of no flare in the solar-wind flux (e.g. Yoshino et al. 2009). 
However, since SWCX and thermal plasma is hardly distin- 
guished with a limited spectral resolution of CCDs (Gupta et 
al. 2009), we model the sum of SWCX and LHB as a single 
unabsorbed thermal plasma. Here, we assumed zero redshift 
and a solar metallicity. We examined the spectral fits with the 
model, by changing the parameters to be free or fixed. We 
also assumed that there was no intracluster and/or intergalac- 
tic emission outside of the virial radii. We discuss possible 
intergalactic emission in subsection 3.5. The resultant param- 
eters of the fits are shown in table 3 and figure 2. As a result, 
the spectra were well-fitted by the sum of Galactic and CXB 
components. We also tested the inferred photon index, F, and 
normalization of the CXB model in Kushino et al. (2002) after 
the point-source subtraction. Our results gave slightly lower in- 
tensity than the Kushino et al. (2002) level, with poor statistics 
above 2 keV. Thus, we employed the parameters of the case 5 
in table 3, derived with fixed spectral shape for the Galactic 
emission and fixed shape and intensity for the CXB compo- 
nent, respectively. We performed spectral fit for the clusters 
in the Sculptor supercluster with this fixed background in the 
following analysis. 

We also examined whether or not an additional apec com- 
ponent (with redshift either zero or supercluster value) gave a 
significant improvement in x 2 against a change of degree of 
freedom Av. In the case of zero-redshift apec model, the im- 
provement was Ax 2 = 1 with Av = 3 over the value of case 4 in 
table 3. This is not a significant change considering the F-test. 
On the other hand, in the case of adding a cluster-redshift apec 
model (z = 0.108), there was no improvement (A% 2 < 1) over 
the case 4 value, for a temperature range between 0. 1-0.8 keV. 

Since the fit is not formally acceptable Of 2 /v = 231/170 in 
case 4), we investigated the x 2 behavior in 2 energy bands, 
0.4-1.0 and 1.0-5.0 keV, separately. As a result, the^- 2 -value 
in 1.0-5.0 keV (the reduced^ 2 ~ 1.5) was much larger than 
that in 0.4-1.0 keV (the reduced^ 2 ~ 1.0). Changing the pa- 
rameters of the CXB component did not improve the fit. The 
spectral data in this high-energy range seem to fluctuate from 
bin to bin by more than the statistical limit. There are some fea- 
tures around the edges of Si-K (1.8 keV) and Au-L (3.4 keV), 
suggesting a problem in the XIS response function. However, 
the model reproduces average feature of the observed spectrum 
quite well, yielding the CXB spectrum consistent with the pre- 
vious Kushino et al. (2002) result, and the fit is good in the en- 
ergy range of our interest (< 1 keV). Therefore, we proceeded 
with the spectral fit using these data. 

We also investigated the effect of possible incorrect calibra- 
tion, i.e., NXB level and OBF contaminations, by artificially 
changing the values by ±10% and comparing the resultant^ 2 
value. While the temperatures of the LHB and MWH did not 
change within < 1 % compared with the values for the case 1 in 
table 3, the normalization and F of the CXB models changed 
by ~10%. We also examined the other uncertainty in the OBF 
contaminant by changing the absorber thickness by +10%. The 
temperatures of LHB and MWH and F and normalization of 
the CXB did not change within < 1%, the normalizations of 
the LHB and MWH changed by ± ~ 10%. A list of * 2 /dof by 
changing the systematic errors in case 4 is presented in table 4. 



We also estimated the fraction of photons entering from out- 
side of the extracted regions using a simulator of the Suzaku 
XRT/XIS system "xissim". The photon contamination can 
cause a significant effect in the case of cluster observations, 
since much bright center is outside of the field and the point- 
spread function of the Suzaku XRT shows an extended tail (see 
e.g. Sato et al. 2007a). If we assume a beta-model of A2811 
emission as described in subsection 2.2 extending beyond the 
virial radius, a half of the simulated photons in the extracted 
region, which is outside of the virial radius, is originated from 
outside of the extracted region. However, this emission has a 
flux of ~ 1/50 of the estimated CXB level below 2 keV. Thus, 
the effect of the stray light is negligible in the spectral fits. 

3.3. Observations close to the Sculptor supercluster 

In order to estimate the Galactic and CXB spectra pre- 
cisely, we looked into other observations which were carried 
out near the Sculptor supercluster. There are two observations: 
NGC 253 offset region, and Q0109-3518. Their observation 
logs are shown in table 5. We analyzed these data with the 
same criteria used for the fore- and background spectral esti- 
mation as described in subsection 3.2. As mentioned in sub- 
section 2.2, the contamination on the XIS optical blocking fil- 
ter has increased after the launch, and the effective area in the 
energy range below 1 keV were smaller in these offset data as 
shown in figure 4. However, because this effect was properly 
included in the ARF, the observations were useful in improving 
the background estimation. 

After subtracting point-like sources as shown in figure 3, 
we carried out spectral fits with the same model, consisting of 
LHB, MWH and CXB, by changing the parameters to be either 
free or fixed. The resultant parameters are shown in table 6 and 
in figure 4. As a result, the spectra of these near-by regions are 
also well-fitted with the background model as shown in subsec- 
tion 3.2. The intensities of the Galactic and CXB components 
are consistent with those for the region outside of the virial ra- 
dius of the clusters in the Sculptor supercluster. 

3.4. Radial Profile oflCM Properties for each cluster 

We extracted spectra from the regions < r < 0.1 riso,0.1 < 
r < 0.3 ri8o,0.3 < r < 0.5 r^o, and 0.5 < r < 1 rigo, centered 
on (RA, Dec) = (0 h 42 m 09. s 3, -28°32'33") for A2811, < r < 
0.4 ri8o,0.4 < r < 0.6 riso, and 0.6 < r < 1 r\w, centered on 
(RA, Dec) = (0 h 38 m 32. s 3, -29°04'44") for A2801, and 0< r < 
0.5 r 180 , 0.5 < r < 1 .0 r 180 centered on (RA, Dec) = (0 h 39 m 38?7, 
-28°53'21") for A2804, respectively. The spectra with BI and 
FI for all regions were fitted simultaneously. Energy ranges 
used for the fit are 0.4-7.1 keV for A2811 and A2801, and 
0.4-5.0 keV for A2804, respectively. In the simultaneous fit, 
the Galactic and CXB components based on the case 4 values 
as described in subsection 3.2 were commonly included for all 
regions. We excluded the narrow energy band around Si K- 
edge (1.825-1. 840 keV) because its response was not modeled 
correctly. The energy range below 0.4 keV was also excluded 
because the C edge (0.284 ke V) seen in the BI spectra could not 
be reproduced well in our data. These are the known problems 
in the XIS response. 

The BI and FI spectra of the ICM in A2801 and A28 1 1 were 
well-represented by a single temperature model as shown in 
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Fig. S. The panels show the observed spectra for the regions of A2801, A2804, and A281 1 which are denoted in the panels, and the data are plotted by 
black and red crosses for BI and FI, respectively. The estimated NXB components are subtracted, and the black and red lines show the best-fit model for 
the BI and FI spectra, respectively. The BI spectra of the ICM component are shown in an orange line. The CXB components are shown in a magenta 
line, and the Galactic emissions are indicated by cyan and blue lines, respectively. For A2804, the additional ICM and LMXB components are shown 
in orange and magenta dashed lines, respectively. The energy range around the Si K-edge (1.825-1.840 keV) is ignored in the spectral fits. The lower 
panels show the fit residuals in units of <x. 



figure 5. The spectra except for the r < 0.3 r^o region of 
A28 1 1 were all fitted with apec model, while these inner re- 
gion spectra of A2811 needed fitting with vapec model. In 
this fit, the metal abundances were linked in the following way, 
Mg=Al, S=Ar=Ca, and Fe=Ni. As for A2804, the spectrum in 
0.5 < r < 1 ri8o was represented by a single temperature model, 
but the data in the inner r < 0.5 r^o region needed an addi- 
tional low mass X-ray binary (LMXB) component represented 
by zbremss with IcT = 10 keV as well as the apec ICM compo- 
nent. The A2804 emission is significantly affected by point- 
like sources due to their lower ICM temperature than those 
in other clusters. Results of the temperature and abundance 
from the spectral fit for individual regions are summarized in 
table 7-8 and figure 6. The typical systematic error on the tem- 
perature of A28 1 1 due to the OBF contamination and change 
of NXB level by +10% are indicated. 

The ICM temperatures of the clusters clearly decrease to- 
wards the virial radii. The temperatures in 0.5-1 nso region 



drop to about 60% of the peak temperature in the three clus- 
ters. This tendency is consistent with the previous Suzaku re- 
sults (Bautz et al. 2009; Reiprich et al. 2009; Hoshino et al. 
2010). While the metal abundances of a elements in the cen- 
tral region of A28 1 1 are about the solar value, the iron abun- 
dance is ~ 0.4 times the solar value (see table 8). The obtained 
temperatures are consistent with the ASCA results (Obayashi 
et al. 1998), but the present results give better accuracy and ra- 
dial distribution of temperature. The metal abundance in the 
central region of A2801 is also 0.4 solar, which is mostly de- 
termined by iron. The metal abundance in the outer region of 
A2801 and the whole region of A2804 had to be fixed at 0.3 so- 
lar because of the poor photon statistics. The iron abundances 
of these clusters in the r < 0.1 rigrj region may be slightly lower 
than those in other clusters (e.g. Tamura et al. 2008; Hoshino 
et al. 2010), while the alpha elements in A2811 show some 
enhanced abundance compared with the iron. 

Kawaharada et al. (2010) reported a higher temperature and 
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Table 7. Summary of the parameters of the fits of A2801 and A2804. Errors are 90% confidence range of statistical errors, and do not include systematic 
errors. The solar abundance ratio of angr was assumed. 



A2801 








Region 


0<r< 0.4 r m 


0.4< r < 0.6 r m 


0.6< r < 1 r m 


kT (keV) 


9 ot+0.15 

-Q45 
43+0.12 

U - H -0.10 


i 43+0.24 


1.08*j3 


Z (solar) 


0.3 (fix) 


0.3 (fix) 


Norm* 


15 9 +L2 
-l.i 


7.01+l^xlO- 1 


6.51^i° x 10~ 2 


X 2 ldof 


211/188 


219/189 


244/189 


A2804 








Region 


0< r < 0.5 r^o 


0.5<r< 1 rigo 




kT (keV) 


1 15 +031 /0 63 +0 23 

laJ -0.13' u -° -0.28 


OQi+0.06 

u -^ 1 -0.07 




Z (solar) 


0.3 (fix) 


0.3 (fix) 




Norm' (1(T 3 ) 


5 37 +102 /l 35 +219 

J ' J -1.87' J -0.79 


2 15 +0 - 27 




X 2 /dof 


176/168 


169/171 





* Normalization of the vapec component scaled with a factor of the selected region comparing to the assumed 
image in "xissimarfgen", Norm = f actor J n^rmdV / [47r(l + z) 2 D^] xl(T 20 cirT 5 arcmin~ 2 , where Da is the 
angular distance to the source. 

t Normalization of the apec component divided by the solid angle, fl u , assumed in the uniform-sky ARF calcula- 
tion (20' radius), Norm = f n e n H dV / (4n(l + z) 2 D 2 )/n u xl0~ 14 cm 4007T arcmin , where Da is the angular 
distance to the source. 



Table 8. Summary of the parameters of the fits of A281 1. Errors are 90% confidence range of statistical errors, and do not include systematic errors. The 
solar abundance ratio of angr was assumed. 



A2811 










Region 


0<r<0.1 riso 


0.1< r < 0.3 r I80 


0.3<r<0.5 riso 


0.5<r< 1 riso 


kT (keV) 


^ to +0.22 

-P.46 


5 17 +al2 

J - L ' -0.12 


4.97!g 


2 63 +0 ' 82 

^■"-'-0.57 


O (solar) 


70+0.60 

-0.56 


f) T T +0.40 

-Q-33 


0.28!°;! 3 


0.3 (fix) 


Ne (solar) 


n qq+0 66 

y -Q.62 


1.08!°* 


4 


I 


Mg, Al (solar) 


0.90!- 


0.80!| 


1 


i 


Si (solar) 


0.85!$ 


0.60!; 


1 


i 


S, Ar, Ca (solar) 


46+0.46 

u ' HU -0.46 


53+ 033 


i 


i 


Fe, Ni (solar) 


34+ 05 


35 +0 04 

u -- ,J -0.04 




i 


Norm* 


648 +3 i 


81.5 +2 - 7 


8.08!°;1 


968 +0169 


X 2 /dof 


194/183 


235/183 


240/188 


215/189 



* Same as the unit of the normalization of A2801 in figure 6 



entropy region in the northeastern edge of A 1689 which was 
connected to an overdense filamentary structure of galaxies 
outside of the cluster. As shown in figure 1 right, the galaxy 
distribution around A28 1 1 shows an extension to the northwest 
direction rather than the southwest direction. We investigated 
the difference of temperature between these directions in the 
radius range of 0.3 < r < 0.5 r^o and 0.5 < r < 1 r^o from 
A281 1. Metal abundance was fixed to 0.3 solar because of the 
poor statistics. As a result, we found no temperature variation 
between the two directions within the statistical error. 

3.5. Upper limit of Oxygen emission lines 

Finally, we fitted the spectra taken outside of the virial radii 
of the clusters. The region is the A28 1 1 offset region, supposed 
to be the filament connecting the clusters. This region is the 
same as the one we used to estimate the fore- and background 
emissions as shown in subsection 3.2. The spectrum in this 
region is well-represented by the Galactic and CXB compo- 



nents without an additional model required as shown in subsec- 
tion 3.2, however we investigate the upper limit of the intensi- 
ties of the additional O VII and O VIII lines at the cluster redshift 
from the WHIM origin. We fitted the data with the following 
model: apec LHB + phabs x (apec MWH + pow CXB + gaussian Q vn + 
gaussian 0vm ). Here, the parameters of the Galactic and CXB 
components followed our previous fit in case 5. The tempera- 
tures of the Galactic components and the photon index of the 
CXB component were fixed, and the normalizations were free 
in the fit. The redshifted central energies of the O VII and O VIII 
gaussian lines were fixed at 518 and 589 eV, respectively, with 
no intrinsic width of the lines assumed. We employed the in- 
crement of Ax 2 — 4 as the measure for the 2<x upper limits of 
the line intensities. This value also corresponds to the 95% 
limit for an additional single parameter (single line intensity) 
in the F test. As a result, we determined the upper limit of the 
O VII and O VIII lines to be 8 . 1 and 5 . 1 x 1 0~ 8 photons cirT 2 s~ 1 
arcmin -2 as the 2<x confidence limits as shown in table 9 and 
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Abell 2811 




Abell 2804 



0.5 I ' 1 ' 1 ' 1 ' 1 ' 1 

0.2 0.4 0.6 0.8 1 

Projected radius (r/r :m ) 

Fig. 6. Radial temperature profiles of A2811, A2801, and 
A2804 against the projected radius in unit of the virial ra- 
dius. Black and Light gray lines show typical systematic 
changes of the best-fit values by varying the thickness of the 
OBF contaminant and NXB levels by ±10%, respectively. 
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Fig. 7. The panels show the resultant fits for the constraint of the inten- 
sities of the O VII and O VIII lines in the filament region between clusters 
of the Sculptor supercluster. The O VII emission line is shown in a green 
line, and the notations of the other lines are same as shown in figure 2. 

figure 7. Even if we allowed the redshift of the O VII and O VIII 
lines to vary to give the peak energies in 490-540 and 560- 
620 eV, respectively, the upper limit did not change within the 
statistical errors. Our resultant values are almost consistent or 
give slightly stronger constraint compared with the results in 
Takei et al. (2007b). This is due to the combination of less 
coupling between the redshifted and the Galactic O lines, to 
slightly lower Galactic foreground intensity and to lower inter- 
stellar absorption. 

In order to constrain the WHIM density based on the O line 
upper limits, we show a wedge diagram of galaxy distribution 
in figure 8. This is expressed in R.A. vs. redshift projected for 
the declination range -32° < Dec. < -26° in the direction of 
the Sculptor supercluster. The galaxy concentration around the 
Sculptor supercluster around z ~ 0.1 1 is obvious and we assume 
~ 25 Mpc for the line-of-sight depth of the WHIM distribution, 
which corresponds to 0.105 <z < 0.1 1 1. As mentioned in Takei 
et al. (2007b), we constrain the density of the WHIM cloud 
assuming the temperature to be T = 2 x 10 6 K. The line in- 
tensity, / < 8.1 x 10~ 8 photons cirT 2 s _1 arcmkr 2 at z=0.1079, 



Table 9. Resultant intensities of O VII and O VIII emission lines with 2<x 
confidence range. 



Ovu 




Center (keV) 


0.518(fixed) 


Sigma 


(fixed) 


I (xl0~ 8 photons cirT 2 s _1 arcmkr 2 ) 


4 2 +3 - 9 


O VIII 




Center (keV) 


0.589(fixed) 


Sigma 


(fixed) 


I (xl0~ 8 photons cirT 2 s arcmkr 2 ) 


+51 



and electron-to-hydrogen number density ratio n e /tin= 1.2 for 
ionized gas give the following condition: 

« H < 1.6 x 10~ 5 cm" 3 (Z/0.1 Z r 1/2 (L/25 MpcT 1/2 . (1) 

The corresponding overdensity 6 = rin/riH is 

5 < 60 (Z/0.1 Z r 1/2 (L/25 Mpc)" 1/2 . (2) 

Compared with the previous Suzaku observations (Takei et 
al. 2007b; Tamura et al. 2008), our upper limits for the O VII 
and Ovm lines imposed slightly tighter constraints on the 
WHIM emission. The present result is consistent with the ex- 
pected values from numerical simulations, indicating a baryon 
over-density p/p C nt ~ 10— 100 by Dolag et al. (2006) for exam- 
ple. We note that the dense part of WHIM is found only in the 
superclusters and the gas is likely to be inhomogeneous with 
a wide range of temperature and density. In the case of the 
bridge region between A222 and A223, Werner et al. (2008) 
showed the temperature to be 0.9 keV. Such a high temper- 
ature (also with high density) yields lower emissivity of O VII 
and O VIII lines. Thus, estimation of temperature and density is 
not straightforward. 

As mentioned in section 1, Fang et al. (2010) reported a 
significant detection of Ovu absorption line in the dkection 
of the "Sculptor Wall" at z ~ 0.03. Zappacosta et al. (2010) 
also searched for the C, O, and Ne absorption lines from the 
"Pisces-Cetus Supercluster" and the " Farther Sculptor Wall" 
in the same line of sight as the "Sculptor Wall". Even though 
the Sculptor supercluster and the above studied filaments are 
not physically connected with each other, we can compare the 
estimated hydrogen density. Our upper limit on the density is 
twice larger than the "Sculptor Wall" value reported by Fang 
et al. (2010). Also, the overdensities derived for the "Pisces- 
Cetus Supercluster" and "Farther Sculptor Wall" are ~ 1/4 
times of our result. This indicates a higher sensitivity of the ab- 
sorption line measurement in the detection of WHIM. In fact, 
Suzaku cannot probably detect the O line emission from the 
"Sculptor wall" because the redshift is too low to separate from 
the Galactic emission lines. Thus, future quantitative studies 
of the WHIM will need X-ray optics with larger effective area 
and instruments with much higher energy resolution such as 
microcalorimeters, which would play a key role in future in- 
vestigation of the cosmic plasma diagnostics. 
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Wedge diagram in R.A. vs. redshift 

Fig. 8. Wedge diagram expressed in R.A. vs. redshift projected on 
-32° < Dec. < -26° . The data are taken from 2dFGRS as mentioned in 
subsection 2. 1 and figure 1 . About 44,000 galaxies exist in the direction 
of the map and the redshift, z < 0.2. The red and green dashed line show 
the direction of our observations and "the Sculptor Wall" at z ~ 0.03 
in Fang et al. (2010), respectively. The red and green arrows also 
show the Sculptor supercluster and "the Sculptor Wall", respectively. 

4. Summary 

Suzaku carried out observations of the Sculptor superclus- 
ter in four fields: A2811, A2811 offset, A2804, and A2801 
regions. After estimation of the foreground and background 
emission, we examined the annular cluster spectra to the virial 
radius. The resultant temperature of the clusters showed sig- 
nificant decrease from the central region toward the virial ra- 
dius, consistent with the previous report for other clusters. Iron 
abundance suggests a nearly uniform distribution. We also 
searched for the oxygen line emission outside of the virial 
radii of the clusters, corresponding to the filament region in the 
Sculptor supercluster. The spectral fits on the filament region 
data did not require an additional component except for the 
Galactic and CXB components. We constrained the intensities 
of redshifted O VII and O VIII lines to be less than 8. 1/5. 1 xl0~ 8 
photons cirT 2 s _1 arcmin~ 2 , respectively, as the 2<x upper lim- 
its. The intensity of O VII line indicates ;?h < 1.6 x 10~ 5 cm -3 
(Z/0.1 Z r 1/2 (L/25 Mpc)- 1/2 , and the over density, 6 < 60 
(Z/0.1 Z r 1/2 (L/25 Mpc)- 1/2 . 
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